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Abstract 

An improved measurement of vi for tt in a broad range of pj and centrality is presented. By 
combining V2 with the Raa, we provide new insights on jet-medium interactions. We show that 
current pQCD energy loss models cannot describe the suppression of the n° as a function of 
the angle with respect to the reaction plane. Our result could help to resolve the factor of 4 
differences in the predicted transport coefficients among these models. Alternatively, it may 
suggest that non-perturbative effects associated with the strongly coupled QGP are important, 
and new theoretical developments are needed to fully understand the jet medium interactions. 

The measurement of n a anisotropy at PHENIX potentially serves three important goals for 
heavy-ion program at RHIC: a) to better constrain the energy loss processes by studying the sup- 
pression of 7T° as function of angle to the reaction plane (RP); b) to investigate the mechanisms for 
anisotropy at intermediate pp, where the transition from a pressure gradient dependent outward- 
push driven by collective flow to a path length dependent attenuation driven by jet quenching 
takes place; c) to understand the non-flow contribution, especially those from jets, in the high 
pj v'2 measurements. These goals are motivated by the current limitation in our understanding 
of the interactions between jets and flowing medium. For example, the extraction of transport 
properties suffer from large theoretical uncertainties. Predictions by various jet quenching frame- 
works differ by a factor of 4, even though they all describe the Raa measurements with identical 
3-D hydro evolution of the underlying medium |l|]. Precision measurements of RP dependent 
suppression can provide severe constraints on these models. 

These goals are facilitated by the large increase in the effective statistics of the 2007 Au+Au 
dataset from PHENIX, which is the basis of current analysis. It has two important improvements 
over the previous PHENIX measurements JUS: 1) various new RP detectors were installed 
covering a broad rj range with improved resolution (RxNPin, RxNPout, MPC at 1.0 < 77 < 1.5, 
1.5 < 77 < 2.8 and 3 < r\ < 4, respectively), and 2) factor of 4 increase in collected statistics 
relative to Jit]. This is equivalent to an effective factor of 14 and 7 increase for tt° measurement 
for RxNP and MPC, respectively. Due to limited space, we shall focus our discussion on goal a) 
and briefly mention goals b) and c) at relevant places. 

The analysis follows the same technique laid out in JH. The main idea is to extract 7r°'s from 
the invariant mass distribution constructed separately in six bins of n° angle relative to the event 
plane (EP), A<p — <p^ - *¥ep, in the interval A<p e [0 - n/2]. The raw V2 is obtained by fitting the 
raw yield with A(l + 2v™ w cos 2A<p), which is then corrected by the corresponding reaction plane 
resolution, cr RP , to give the signal V2 — v™ v I<Jrp- Note that the EP is obtained via standard event 
plane method and cr RP is measured via the sub-event method by correlating the north and south 
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RP detectors (which are symmetric). The differential nuclear modification factor Raa(A<Pi, pj) is 
calculated from the published angle-averaged Raa(Pt) II] as 



Raa(M>uPt) = Raa(pt) 



N(Atpi, p T ) 1 + 2v2 cos 2A<fo 
< N(A(f>, p T )>l+ 2v r 2 aw cos 2A0 



(1) 



One advantage of the broad r\ coverage for the RP detectors is that one can evaluate the 
potential non-flow effects. Figure [T| shows the V2 vs N par , in a low pj region and a high pj 
region, measured by different RP detectors. There is a clear r\ dependence of V2 values at high 
pj, especially in the peripheral collisions. This is not the case at low pj. This suggests that most 
of the non-flow effects are induced by jets, and they only become important at high p j and for RP 
detectors reside closer to central arm where the 7r°s are detected. Even though RxNP detectors 
have the best resolution, they may suffer jet bias. Instead, we use MPC for EP determination. 
The MPC sits at the same r\ range as BBC detectors (used in previous analysis); however it has 
40% better resolution compare to the BBC since it measures both neutral and charged particles, 
and has higher granularity. 
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Figure 1: v 2 vs N part 
at a low pj and a high 
pj regions using differ- 
ent reaction plane de- 
tectors. 
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Figure [2] shows V2 values measured in six centrality classes, overlaid with the published 
results from 2004 Au+Au run yfl. The two measurements are consistent, but new data improve 
significantly on both statistical errors and pj reach. The V2 at pj > 6 GeV/c, for all centralities, 
remains significantly above zero and is constant with pj. 

If the particle production at high pj (> 6 GeV/c) is dominated by fragmentation of the jets 
that survive the medium, as most energy loss models predict, then it would lead to an anisotropy 
of the suppression which can be compared with our data. Figure[3]shows the high pj Raa for jt° 
measured in-plane (0 < A<p < n/12) and out-of -plane (5tt/12 < A<p < 6n/l2) directions, derived 
according to Eq. Q] A clear split of the suppression levels can be seen between the two for 
three centrality bins shown (20-30%, 30-40% and 40-50%). Calculations from three pQCD jet 
quenching models (abbreviated as AMY, HT, ASW) from [1] are compared with our data. These 
calculations are carried out with identical initial conditions, medium evolution and fragmentation 
functions, and their suppression levels have been tuned to reproduce the inclusive Raa data in 
central Au+Au collisions as well as the 20-30% bin. However, as Figure [3^ shows, all three 
models under-predict the difference between in-plane and out-of-plane Raa- The AMY model 
describes the out-of-plane Raa, but not the in-plane; HT model is the opposite; ASW model 
does a reasonable job for both in- and out-of plane directions at high pj, but misses the low pj 
region. However, we are not advocating that one model is better than the other, because changing 
other control parameters, like initial conditions, medium evolution or fragmentation functions, 
may change this comparison. Our principal conclusion is that by utilizing both observables, one 
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Figure 2: 7r° V2 vs in six centrality bins. The solid points are the PHENIX preliminary results in Run 7 measured 
with the MPC detector. The open points are the published results from Q|. 



should be able to narrow down the uncertainties on initial condition and medium evolution and 
other effects that were ignored so far in the calculation, e.g. virtuality difference between the 
in- and out-of plane. Our data may help to resolve the discrepancies, e.g. q, between the model 
predictions. 

On the other hand, the failure of these models to reproduce the large anisotropy at high pj 
may imply that the pQCD treatment of the energy loss process as sequential radiation, being pro- 
portional to the local color change density is not sufficient. In the presence of the strongly cou- 
pled medium (sQGP), both the path length dependence and the color charge dependence could 
be modified. For example, calculation based on ADS/CFT technique suggests that AE cc L? ^ 
and q cc yjasYMN c [6] instead of AE cc L 2 and q cc a s N 2 for pQCD. In a separate calculation, 
Kharzeev QVJ] estimates dE/dx cc E 2 in strong coupling limit, much stronger than the logarithmic 
dependence expected from pQCD. These non-linear dependences could be the reason for the 
large anisotropy. Figure [4] shows the comparison of v>2 with several toy model predictions with 
different medium density dependence. The model from [8] introduces strong non-linear depen- 
dence by assuming that jet quenching happens mostly close to phase transition boundary. The 
model from [9] does so by suppressing the energy loss at high energy density region by using a 
plasma formation time argument. Both are able to qualitatively describe centrality dependence. 
However, much need to be done to generalize these toy models into rigorous theoretical calcula- 
tions. 

In summary, we presented high statistics measurement of n° V2 and Raa(A<P, Pt) up to 13 
GeV/c in broad centrality ranges of Au+Au collisions at y saw =200 GeV. Current energy loss 
models were not able to describe Raa(A<P, Pt) at pj- > 6 GeV/c. This discrepancy implies that the 
current pQCD models need further tuning, in which case our measurement can help to resolve the 
differences in these models. It may also suggest that the non-perturbative effects associated with 
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Figure 3: Data points: solid(open) is the in-plane(out-of-plane) n° Raa(Pt) for three centrality bins measured with the 
MPC detector. Lines represent calculations of AMY, HT and ASW models from [1], respectively. 



Figure 4: Data points: V2 val- 
ues at pr >6GeV/c measured with 
the MPC detector. Open sym- 
bols: pQCD model calculations 
taken from [II]. Lines: geometric 
model calculations with different 
assumptions on path length depen- 
dence OSS. 
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the strongly coupled QGP are important, and new tools are required to understand the energy loss 
processes. Further detailed study of the V2 and R&a in the intermediate pj range of 2-6 GeV/c 
may shed light on the role of perturbative and non-perturbative effects. 
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